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The influence of the finite temperature on the antiferromagnetic ( AFM) spin ordering in symmetric 
nuclear matter with the effective Gogny interaction is studied within the framework of a Fermi liquid 
formalism. It is shown that the AFM spin polarization parameter of partially polarized nuclear 
matter for low enough temperatures increases with temperature. The entropy of the AFM spin 
state for some temperature range is larger than the entropy of the normal state. Nerveless, the free 
energy of the AFM spin state is always less than the free energy of the normal state and, hence, the 
AFM spin polarized state is preferable for all temperatures below the critical temperature. 

PACS numbers: 21.65.+f; 75.25.+z; 71.10.Ay 



INTRODUCTION 

The spontaneous appearance of spin polarized states 
in nuclear matter is a topic of a great current interest 
due to its relevance in astrophysics. In particular, the 
effects of spin correlations in the medium strongly influ- 
ence the neutrino cross section and neutrino mean free 
path. Therefore, depending on whether nuclear matter is 
spin polarized or not, drastically different scenarios of su- 
pernova explosion and cooling of neutron stars can be re- 
alized. Another aspect relates to pulsars, which are con- 
sidered to be rapidly rotating neutron stars, surrounded 
by strong magnetic field. One of the hypotheses to ex- 
plain such a strong magnetic field of a pulsar is that it 
can be produced by a spontaneous ordering of spins in 
the dense stellar core. 

If a spin polarized state appears, nucleons with spin up 
and spin down lie on different Fermi surfaces. Spin-spin 
and spin-isospin correlations, leading to the formation of 
a spin polarized state, essentially depend on the overlap 
between Fermi surfaces. It is controlled by the number 
densities of nucleons of different species and spin polar- 
ization parameter, which, in turn, should be determined 
self-consistently. 

From the general point of view, the problem of finding 
the phase diagram of a many-particle Fermi system when 
fermions lie on different Fermi surfaces is encountered in 
many physical cases. For example, thermodynamic prop- 
erties of a neutron-proton condensate in asymmet- 
ric nuclear matter are governed by pairing correlations 
between neutrons and protons, occupying two different 
Fermi spheres. Analogous situation appears in high den- 
sity QCD, when the quark system is unstable against 
the formation of the (q q) color superconducting conden- 
sate 0-E3- Besides, one can mention the spin singlet 
pairing in a superconducting metal in the presence of 
magnetic impurities 0, or the B-phase of the su- 
perftuid 3 He (Balian-Werthamer phase) in a magnetic 
field [13j, when the spin degeneracy is loosened due to 
the Pauli paramagnetism. Another example is the ap- 
pearance of BCS pairing in ultracold trapped quantum 



gases El El [13. 

The possibility of a phase transition of normal neu- 
tron and nuclear matter to the ferro mag netic (FM) spin 
state was studied by many authors |l7 | -|26 | . predicting 
the ferromagnetic transition at g w (2-4) go for differ- 
ent parametrizations of Skyrme forces (go = 0.16 fm -3 is 
the nuclear matter saturation density). Competition be- 
tween FM and AFM spin ordering in symmetric nuclear 
matter with the Skyrme effective interaction was studied 
in Ref . [27| , where it was clarified that the FM spin state 
is thermodynamically preferable to the AFM one for all 
relevant densities. However, strongly asymmetric nuclear 
matter with Skyrme forces undergoes a phase transition 
to a state with oppositely directed spins of neutrons and 
protons |28j . The same conclusion in favour of antipar- 
allel ordering of neutron and proton spins in symmetric 
nuclear matter was confirmed also in Ref. [29] for the 
Gogny effective interaction, where it was shown that the 
AFM spin state appears at g « 3.8go- 

For the models with realistic nucleon-nucleon (NN) in- 
teraction, the ferromagnetic phase transition seems to be 
suppressed up to densities well above go [3(|-[32- In 
particular, no evidence of ferromagnetic instability has 
been found in recent studies of neutron matter |33j and 
asymmetric nuclear matter ,34] within the Brueckner- 
Hartree-Fock approximation with realistic Nijmegen II, 
Reid93, and Nijmegen NSC97e NN interactions. The 
same conclusion was obtained in Ref. |35|. where the 
magnetic susceptibility of neutron matter was calculated 
with the use of the Argonne vis two-body potential and 
Urbana IX three-body potential. 

Here we continue the study of spin polarized states in 
nuclear matter, using as a potential of NN interaction 
the effective Gogny forces |H E3 and assuming that the 
AFM spin ordering is realized as a ground state of nuclear 
matter at zero temperature |29j | . The main emphasis will 
be laid on determining the finite temperature behavior 
of the AFM spin polarization. The first goal of the study 
is to show that at low enough temperatures thermal fluc- 
tuations promote the AFM spin polarization of nuclear 
matter, when a system of nucleons can be treated as a 
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multicomponent Fermi liquid |38l l39l |40j. The second 
goal is to provide a fully self-consistent calculation of 
the basic thermodynamic functions of antiferromagneti- 
cally ordered nuclear matter at finite temperatures with 
a modern effective finite range NN interaction. In spite 
of that the entropy of the AFM spin state can be larger 
than the entropy of the normal state, the total balance 
of free energies lies with the AFM spin state for all tem- 
peratures below the critical temperature. 



BASIC EQUATIONS 

The normal states of nuclear matter are described by 
the normal distribution function of nucleons / KlK2 = 
Tr ga^ a Kl , where n = (p, a, r), p is the momentum, 
<t(t) is the projection of spin (isospin) on the third axis, 
and g is the density matrix of the system. Bearing in 
mind to consider the possibility of FM and AFM phase 
transitions, the normal distribution function / and the 
nucleon single particle energy e can be expanded in the 
Pauli matrices Oi and in spin and isospin spaces 

/(P) = /00(P)^0T + i3o(p)cr 3 To (1) 
+ /03(P)OOT3 + /33(P)CT3T3- 

e(p) = £oo(p)ctot + £3o(p)o" 3 to (2) 

+ £03(p)o- T 3 + £33(P)C-3T3- 

Expressions for the distribution functions /oo, /30, /o3> /33 
in terms of the quantities e read [271 128| 

/oo = -An{u n {) + "-(^pt) + + n(uj pl )}, 
ho = T{n(wn-r) + n{u P i) ~ n(u nl ) - n(w p ;)}, (3) 
/os = T{"(wnr) - n{u P i) + n(u nl ) - n(wpi)}, 
h-3 = -An{u n -\) - n(u p] ) - n(uj nl ) + n(uj pl )}. 
Here n(ui) = {exp(w/T) + l}" 1 and 

w nT = CoO + £30 + £03 + £33i 

w pT = Coo + £30 — £03 — £33 , (4) 
w n| = Coo — C30 + ?03 — ?33j 
Vpl = COO - £30 - £03 -I- £33 , 

where 

Coo = £00 — Moo, £30 = e 3 o, 
£03 = £03 — M03 > £33 = £33 ! 

Mti ~r Mp Mti Mp 

Moo = 2 ' = 2 ' 

/i„, fi p being the chemical potentials of neutrons and pro- 
tons. The branches u> n ^,uj n ^ of the quasiparticle spec- 
trum correspond to neutrons with spin up and spin down, 



and the branches u) p ^ , u P i correspond to protons with 
spin up and spin down. 

The distribution functions / should satisfy the normal- 
ization conditions 



V 



^2ho(p) = Q, 



y ^2 /o3(p) = Qn ~ Qp = ag, 



V 



X! /so(p) = Q\ - Qi = Ag TT , 



(5) 
(6) 



y /33(p) = {Qw\ + Qpi) - {Qni + £>pt) = A en- ( 8 ) 



Here a is the isospin asymmetry parameter, g>„j, g n ± and 
0P? > Spl are the neutron and proton number densities 
with spin up and spin down, respectively; g-\ = g n -\ + g p -\ 
and gi = g n ± + g P i are the nucleon densities with spin 
up and spin down. The quantities Agjj and Agyplay 
the roles of FM and AFM spin order parameters [28| . 

The self-consistent equations for the components of 
the single particle energy have the form [27I l2g| 

£00 (p) = £o(p) + £00 (p) - Moo, 60 (p) = £30 (p), (9) 
£03 (P) = £03 (P) ~ M03, £33 (P) = £33 (P)- 

Here £o(p) is the free single particle spectrum, and 
£00 j £30 j £03 j £33 are the FL corrections to the free single 
particle spectrum, related to the normal FL amplitudes 
C/o(k), C/ 3 (k) by formulas 



p q 



(10) 



£oo(p) = ^^C/o(k)/oo(q), k 
q 

Mp) = -^X>i(k)/ao(q), 
q 

eo3(p) = ^E^2(k)/o3(q), 
q 

£33(p) = i^%(k)/ 33 (q). 



To obtain numerical results, we use the effective Gogny 
interaction D1S 37] . Expressions for the normal FL am- 
plitudes in terms of Gogny force parameters were writ- 
ten in Ref. [2^. Thus, with account of expressions © for 
the distribution functions /, we obtain the self-consistent 
equations © , (|10|l for the components of the single parti- 
cle energy £ o(p), 60 (p), £03 (p), £33 (p), which should be 
solved jointly with the normalization conditions JSJ-JHJ), 
determining the chemical potentials Moo,Mo3, FM and 
AFM spin order parameters Agfj, Agfj_. 

To examine the thermodynamic stability of different 
solutions of self-consistent equations, it is necessary to 
compare the corresponding free energies F = E — TS, 
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FIG. 1: AFM spin polarization parameter as a function of 
temperature at different densities for the D1S Gogny force. 

where the energy functional E is characterized, in the 
general case, by four FL amplitudes Uq,...,U% and the 
entropy reads 

S = J! { n i UJ ra) Inn (uJra) 

p T=n, p cr=T, | 

+ n(uj Ta ) lnn(w rcr )}, n(u) = 1 - n(ui). 

PHASE TRANSITIONS AT FINITE 
TEMPERATURE 

Further we will consider symmetric nuclear matter 
(Qn = Qp)- It was shown in Ref. |29| that in symmet- 
ric nuclear matter with D1S Gogny interaction only the 
AFM spin ordering is realized at zero temperature, but 
the self-consistent equations have no solutions at all cor- 
responding to the FM spin ordering. Our aim here is 
to study the temperature behavior of the AFM spin po- 
larization in the whole temperature domain below the 
critical temperature, T <T C . 

In the AFM spin state of symmetric nuclear matter 
Qnl = QpliQnl = f?p|' neutrons with spin up and pro- 
tons with spin down fill the Fermi surface of radius ki 
and neutrons with spin down and protons with spin up 
occupy the Fermi surface of radius k\ , satisfying at zero 
temperature the equations 

-^[k\ - k\) = Ag n , JL(fc3 + fc3) = e . 

Now we present the results of the numerical solution 
of the self-consistent equations with the D1S Gogny ef- 
fective force. In Fig. 1 it is shown the dependence of the 
AFM spin polarization parameter Ag^/g as a function 



FIG. 2: The free energy per nucleon, measured from its value 
in the normal state, for the AFM spin state as a function of 
temperature at different densities for the D1S Gogny force. 

of temperature at different fixed densities. The interest- 
ing feature is that, if at T = we have partially AFM 
polarized state (Aqu/q < 1), then under increase of tem- 
perature within some temperature interval the AFM spin 
polarization parameter increases as well. This behavior 
is in contrast with the intuitive supposition that thermal 
fluctuations act as a destroying factor on spin ordering. 
Oppositely, for not too large temperatures thermal fluc- 
tuations in nuclear matter promote the AFM spin order- 
ing. The reason for such a behavior is that thermal fluc- 
tuations smear Fermi surfaces of nucleons, leading, thus, 
to increasing overlap between Fermi surfaces. Since in- 
teraction between free nucleons is most strong between 
a neutron and a proton in the spin triplet state, then, 
mainly, due to this interaction, modified by the medium, 
some of neutrons with spin down and protons with spin 
up undergo spin flip transitions from the inner Fermi sur- 
face to the outer one. Thus, initial increase of AFM 
spin polarization with temperature is a result of influ- 
ence of thermal effects and medium correlations. Under 
further increasing temperature thermal fluctuations sup- 
press AFM spin ordering, until it completely vanishes. 

In Fig. [3 the difference between the free energies per 
nucleon of the spin ordered and normal states is shown 
as a function of temperature for different fixed densities. 
One can see that, first, as a result of the initial increase 
of AFM spin polarization, the free energy of the AFM 
spin state decreases with temperature and after that the 
difference between the free energies of the AFM and nor- 
mal states becomes smaller, until it vanishes at critical 
temperature, dependent on density. 

Unexpected moment appears if we consider separately 
the temperature behavior of the entropy of the AFM 
state. In Fig. [3| the difference between the densities of 
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FIG. 3: Same as in Fig. 2, but for the density of entropy, 
measured from its value in the normal state. 

the entropies of the AFM and normal states is shown as 
a function of temperature. One can see that for low tem- 
peratures the entropy of the AFM state is larger than 
the entropy of the normal state. It looks like the AFM 
state at low finite temperatures is less ordered than the 
normal state. Under further increasing temperature the 
difference of the entropies changes its sign and becomes 
negative, that corresponds to the intuitively expected be- 
havior. In spite of that the entropy of the AFM state may 
be larger or less than the entropy of the normal state, 
the difference of the free energies preserves its sign for all 
temperatures below T c . Note that analogous anomalous 
behavior of the entropy with temperature was observed 
also in superfluid asymmetric nuclear matter || . In that 
case for low enough temperatures the entropy of the su- 
perfluid state is larger than that for the normal state. 
Nerveless, as in our case, the total balance of the free 
energies is preserved in favour of superfluid state for all 
temperatures below the critical temperature. It is worthy 
to note that for the AFM spin state anomalous behavior 
of the entropy is observed already in symmetric nuclear 
matter while in the superfluid case it is observed only 
at finite isospin asymmetry. The difference is that in the 
spin ordered state the separation of Fermi surfaces is con- 
trolled by the AFM spin polarization parameter which, 
in turn, is not an independent quantity, but should be 
found self-consistently. 

Note that the stability of the equation of state of nu- 
clear matter with the Skyrme effective interaction in 
terms of Landau parameters was examined in Ref. |4l| . 
where the stability conditions were formulated as the in- 
equalities for the Skyrme force parameters. This study is 
based on the approximation of the effective mass, when 
the quadratic terms on momentum in the Skyrme interac- 
tion are incorporated in the single particle spectrum. The 



approximation of the effective mass, being independent 
of temperature as in Ref. |4l| , is a strong simplifying as- 
sumption and cannot explain the change in the sign of the 
difference between entropies of polarized and unpolarizcd 
states at certain temperature, as seen from Fig. 3. In the 
general case of a finite range interaction, like Gogny force 
in our case or Paris NN potential in Ref. [jj, the single 
particle spectrum is to be determined self-consistently by 
solving the corresponding integral equations (in our case, 
Eqs. Thus, the anomalous behavior of the entropy 

with temperature should be associated with the compli- 
cated renormalization of the free single particle spectrum 
in a strongly interacting nucleon medium. 

Note that, since in the present study only symmet- 
ric nuclear matter has been considered, the obtained re- 
sults cannot be directly extrapolated to proto-neutron 
stars, whose core represents essentially asymmetric nu- 
clear matter. For strongly isospin asymmetric system, 
the appropriate choice of the Gogny force is the use of 
the DIP parametrization ^3|> giving the correct behavior 
of the energy per nucleon at high densities. 

In summary, we have considered AFM polarized states 
in symmetric nuclear matter with the effective Gogny 
interaction at finite temperatures. It has been shown 
that the AFM spin polarization initially increases with 
temperature as a result of smearing Fermi surfaces of nu- 
cleons due to thermal fluctuations and spin and isospin 
dependent correlations in the medium. While the dif- 
ference of the entropies of the AFM and normal states 
anomalously changes its sign at certain temperature, the 
total balance of the free energies lies with the AFM spin 
state for all temperatures below the critical temperature. 



[1] T. Aim, G. Ropke, and M. Schmidt, Z. Phys. A 337, 355 
(1990). 

[2] B.E. Vonderfecht, C.C. Gearhart, W.H. Dickhoff, A. 

Polls, and A. Ramos, Phys. Lett. 253B, 1 (1991). 
[3] M. Baldo, I. Bombaci, and U. Lombardo, Phys. Lett. 

283B, 8 (1992). 
[4] T. Aim, G. Ropke, A. Sedrakian, and F. Weber, Nucl. 

Phys. A406, 491 (1996). 
[5] A. Sedrakian, and U. Lombardo, Phys. Rev. Lett. 84, 

602 (2000). 

[6] A. A. Isayev, Phys. Rev. C 65, 031302 (2002). 

[7] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. 

422B, 247 (1998). 
[8] R. Rapp, T. Schafer, E.V. Shuryak, and M. Velkovsky, 

Phys. Rev. Lett. 81, 53 (1998). 
[9] D. Blaschke and CD. Roberts, Nucl. Phys. A642, 197 

(1998) . 

[10] R. Pisarski, and D. Rischke, Phys. Rev. D 60, 094013 

(1999) . 

[11] A.I. Larkin, and Yu. N. Ovchinnikov, Zh. Eksp. Teor. 

Fiz. 47, 1136 (1964) [Sov. Phys. JETP 20, 762 (1965)]. 
[12] P. Fulde and R. A. Ferrell, Phys. Rev. 135, 550 (1964). 
[13] A. Leggett, Rev. Mod. Phys. Rev. 41, 331 (1975). 



■5 



[14] C.A. Regal, M. Greiner, and D.S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 
[15] J. Kinast, S.L. Hemmer, M.E. Gehm, A. Turlapov, and 

J.E. Thomas, Phys. Rev. Lett. 92, 150402 (2004). 
[16] M.W. Zwierlein, C.A. Stan, C.H. Schunck, S.M.F. Rau- 

pach, A.J. Kerman, and W. Ketterle, Phys. Rev. Lett. 

92, 120403 (2004). 
[17] M.J. Rice, Phys. Lett. 29A, 637 (1969). 
[18] S.D. Silverstein, Phys. Rev. Lett. 23, 139 (1969). 
[19] E. 0stgaard, Nucl. Phys. A154, 202 (1970). 
[20] A. Viduarre, J. Navarro, and J. Bernabeu, Astron. As- 

trophys. 135, 361 (1984). 
[21] S. Reddy, M. Prakash, J.M. Lattimer, and J. A. Pons, 

Phys. Rev. C 59, 2888 (1999). 
[22] A.I. Akhiezer, N.V. Laskin, and S.V. Peletminsky, Phys. 

Lett. 383B, 444 (1996); JETP 82, 1066 (1996). 
[23] S. Marcos, R. Niembro, M.L. Quelle, and J. Navarro, 

Phys. Lett. 271B, 277 (1991). 
[24] T. Maruyama and T. Tatsumi, Nucl. Phys. A693, 710 

(2001). 

[25] M. Kutschera, and W. Wojcik, Phys. Lett. 223B, 11 
(1989). 

[26] A. Beraudo, A. De Pace, M. Martin i, and A . Molinari, 
Annals Phys. 311, 81 (2004); Arxiv: nucl-th/0409039 

[27] A. A. Isayev, JETP Letters 77, 251 (2003). 

[28] A. A. Isayev, and J. Yang, Phys. Rev. C 69, 025801 
(2004). 



[29] A. A. Isayev, and J. Yang, Phys. Rev. C 70, 064310 
(2004). 

[30] V.R. Pandharipande, V.K. Garde, and J.K. Srivastava, 

Phys. Lett. 38B, 485 (1972). 
[31] S.O. Backmann and C.G. Kallman, Phys. Lett. 43B, 263 

(1973). 

[32] P. Haensel, Phys. Rev. C 11, 1822 (1975). 
[33] I. Vidaha, A. Polls, and A. Ramos, Phys. Rev. C 65, 
035804 (2002). 

[34] I. Vidana, and I. Bombaci, Phys. Rev. C 66, 045801 
(2002). 

[35] S. Fantoni, A. Sarsa, and E. Schmidt, Phys. Rev. Lett. 

87, 181101 (2001). 
[36] J. Decharge and D. Gogny, Phys. Rev. C 21, 1568 (1980). 
[37] J.F. Berger, M. Girod and D. Gogny, Comp. Phys. 

Comm. 63, 365 (1991). 
[38] A.I. Akhiezer, V.V. Krasil'nikov, S.V. Peletminsky, and 

A. A. Yatsenko, Phys. Rep. 245, 1 (1994). 
[39] A.I. Akhiezer, A.A. Isayev, S.V. Peletminsky, A.P. 

Rekalo, and A.A. Yatsenko, JETP 85, 1 (1997). 
[40] A.I. Akhiezer, A.A. Isayev, S.V. Peletminsky, and A.A. 

Yatsenko, Phys. Rev. C 63, 021304(R) (2001). 
[41] J. Margueron, J. Navarro, and N. V. Giai, Phys. Rev. C 

66, 014303 (2002). 
[42] M. Farine, D. Von-Eiff, P. Schuck, J.F. Berger, J. 

Decharge, and M. Girod, J. Phys. G 25, 863 (1999). 



